Boson peak in an harmonic scalar model 
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We study, analytically and numerically, an ofT-lattice model of scalar harmonic vibrations for 
structural glasses. The model has a Boson-Peak which we argue can be considered as a prototype 
for materials with a Boson Peak frequency that decreases with lowering temperature. The density 
evolution of the Boson Peak in silica is qualitatively reproduced. The dispersion relation is linear 
at the Boson Peak frequency, in agreement with experiments. In our model the Boson Peak is a 
precursor of a nearby mechanical instability. The Boson Peak is built up by the hybridization of the 
sound waves with extended, but non propagating (in the sense of a fiat dispersion relation), modes. 

PACS numbers: PACS 61.43.Fs, 63.50.+X 



The vibrational dynamics of glasses at the frequencies 
that control the specific heat in the 1 K-lOO K range is 
a long-standing physical problem. At the lowest tem- 
peratures, an anomalous linear dependence on T can be 
reproduced by the two-level model [Q . At slightly higher 
temperatures, another deviation from the Debye law 
is found: the specific heat divided by shows a peak 
near lOK This peak points to the existence of an ex- 
cess of vibrational states at frequencies co 1 THz, which 
shows as a peak in the plot of g(uj)luP' {g{uj) being the 
vibrational density of states, as obtained by Raman or in- 
elastic neutron scattering). Since the scattering intensity 
at the peak scales in temperature with Bose statistics, 
the peak has become known as the Boson peak (BP). The 
physical origin of the BP is a matter of a lively debate 
within the experimental and theoretical glass community. 

The recent experimental efforts to elucidate the nature 
of the high frequency excitations in glasses have produced 
a wealth of inelastic scattering data (X-ray, Raman and 
neutrons). These experiments have found a BP on almost 
all studied glassformers, although of an intensity that 
can vary substantially among different materials. Two 
features of the BP should be remarked. First, the peak 
frequency wbp is several times smaller than any natu- 
ral frequency scale, like the Debye frequency or the band 
edge. Moreover, where good data for the dispersion rela- 
tion Lo^p) (determined from the the Brillouin peak in the 
dynamic structure factor, which probes only longitudi- 
nal modes) are available, it has been checked that uj{p) is 
still a linear function of the momentum p at cjbp • This is 
the case for glycerol, LiCl:6H20 Cao.4Ko.6(N03)i.4, 
B2O3 [D, silica Q and polybutadiene |^ (the interpre- 
tation of scattering experiments in silica has been con- 
troversial 0). Second, although less experimental data 
are available, it seems to be the rule that ubp shifts 
to lower frequencies on heating, while the BP intensity 



grows. This is the case in polybutadiene, polystyrene 
LiCl Q and B2O3 Silica, however, shows the oppo- 
site behavior |llj. Yet, in silica the BP evolution upon 
increasing the density has been studied experimentally 
and in simulations [|l3j. In close agreement with the 
results presented here, it was found that ubp shifts to 
larger frequencies and the BP looses intensity when the 
density grows. The BP has also been identified in the 
dynamic structure factor at high exchanged momentum, 
both experimentally fl^ and in simulations ||l5|] , which 
reflects the fact that the dynamic structure factor tends 
to the density of state (DOS) at high momentum ||l^, 

Perhaps the only widely-agreed upon statement about 
the BP is that it, as well as the whole dynamic struc- 
ture factor in this frequency region, can be understood 
invoking only harmonic vibrations M H, 



|23| , 24 1 (see however the work on soft-potentials [25| for 
a dissenting view). Otherwise, the BP has been vari- 
ously interpreted as arising from mixing of longitudinal 
and transversal modes p^, pO| , pGf , hybridization of 
optical and acoustic modes 1 19[, a combination of the 
level repulsion due to disorder and van Hove singulari- 
ties , an associated mechanical instability [|l^ |l] (as 
suggested by the temperature dependence of ljbp), the 
presence of a loffe-Regel crossover at wbp (2^, ^ , or from 
the scattering of sound-waves with localized anharmonic 
vibrations jll , 25 . The BP has also been identified with 
the anomalous oscilation peak found in Mode Coupling 
Theory as adapted to the non-ergodic phase 

In this Letter we discuss analitically and numerically 
an off-lattice model of scalar harmonic vibrations, which 
is a problem on random-matrix theory |p9| , showing that 
it has a Boson peak with the above-mentioned exper- 
imentally identified characteristics. Thus it is unnec- 
essary to invoke anharmonicity or transverse modes to 
explain the BP (although they are certainly present in 
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real glasses), which instead, as we will show, arises from 
the hybridization of extended but non-propagating modes 
with sound waves, induced by a mechanical instability. 
Our model |]l6[ |3^ consists in particles oscillating 
harmonically around (disordered) equilibrium positions, 
along a fixed direction u (i.e. we neglect all transverse 
modes): Xi{t) = x"'^ + U(pi{t), 1,...,N. The poten- 
tial energy is V = - - ^j? = 
"^fj fiMijipj, where / is the second derivative of the pair 
potential and the dynamical matrix matrix (Hessian) is 

N 

M,, = s,, f^^T - - /(^r - ^f)- (1) 
fc=i 

The model is completely defined when / and the distri- 
bution of the x'^'^ are chosen. We take the latter uniform 
in the whole volume, while / is taken to be regular at 
short distances. As discussed in Ref. this is a good 
first approximation to the correlations of the x'^'^. We 
shall consider the family of functions 

/„(r) = (l-arVa^)e-^^/(2-^), (2) 

where < a < 0.2 (the upper bound has to be im- 
posed in order to guarantee a positive sound velocity). 
When a = (Gaussian case), the hessian is strictly pos- 
itive. When a > 0, we have a stable elastic solid at high 
particle-number density, p, while at low enough p, typi- 
cal interparticle distances will be large enough to allow 
negative eigenvalues (imaginary frequencies). Therefore 
the density controls the appearence of a mechanical in- 
stability in this model: there is a phase-transition from 
a region where the short time dynamics is ruled by a 
positive hessian (high p) to the one where an extensive 
number of negative eigenvalues is found (low p). The 
counterpart in real glasses of that transition is the mode- 
coupling transition [|3^ which is rather a crossover, but 
it has been recently argued to correpond to a well 
defined phase-transition on the topological properties of 
the potential energy landscape. It must also be said that 
in disordered lattices models ^ , an instability anal- 
ogous to ours is found if the disorder is strong enough. 
Nevertheless, the bump in g{uj)/uP' found on the lattice 
is intrinsically tied to a characteristic frequency (the van- 
Hove singularity, see ^ and below). 

We first compute the DOS of our model numerically 
for a = 0.1, using the method of moments with a 
box of side L — 128 u (more than 5 x 10^ particles), 
which allows to reconstruct the spectrum up to very low 
frequencies. The critical p is difficult to locate, but at 
p = 0.05cr~^ imaginary frequencies are clearly found. In 
Fig. |l| (top) we show the DOS divided by uP' at several 
densities. We have supressed the smallest frequencies, 
which are unreliable due to finite-size effects, and used 
the reduced form |l^ p^ ^ giio) j liP' vs. ioj ^ in order to 
compare different densities on the same plot. A peak 
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FIG. 1: Top: BP raising in our model with decreasing p; 
the horizontal lines indicate the Debye value. Bottom: nu- 
merical dispersion relation from the dynamic structure factor 
{g[uj)/Lu'^ rotated 90 degrees is also ploted). Lj{p) is linear at 
ll>bv and gets flat at the high-frequency bump in g{ijj)/ijp. 

in g{ijj)/uP' arises on approaching the instability (i.e. on 
lowering the density). As p is reduced, the peak grows 
(relative to the Debye value, also plotted) and moves to 
lower u). The peak lies always in the linear region of the 
dispersion relation (Fig. |l|, bottom). Since decreasing p 
plays the role of increasing temperature in our model, the 
experimental features of the BP are reproduced, without 
transverse or optical modes or anharmonicity, as claimed. 

Let us discuss analytically the BP found in our 
model. Within our random-matrix approach |l6| 
the DOS is obtained from the resolvent G(p, z) = 
{I IN) exp[ip • (xf - xl'Mz - M)-i],fe, where the 
overline stands for the average over a;°^. Defining G{z) = 
G{p = (X),z), the DOS is g{uj) = -^lm.g{u^ + iO+). 
Expanding the resolvent in 1/p and resumming an 
infinite subset of diagrams |17|| , a non-linear equation is 
found for G{p, z) that in the p — > oo limit yields 

^ = - - /(O) - Ag{z) [^p{q)Giq, z), (3) 

where A = {2tt)~^ ^ p {q) d'^ q and f{q) is the Fourier 
transform of f{r). With this equation, one needs to 
know the resolvent at all q to obtain the DOS, due to 
the last term in the r.h.s. This can be done by solving 
numerically the self-consistent equation of ref. Jl^l • Here 
we perform an approximate analysis, which is more illu- 
minating. The crudest approximation is to neglect this 
term, in which case Eq. ^ is quadratic in Q, and one eas- 
ily finds a semicircular DOS, with center at lo'^ — pf{0) 
and radius 2^/ pA. Indeed, when p oo, the spectrum 
is made of plane waves, with dispersion relation uP'{p) = 
p{f{0)~f{p)), jl^, im (a continous elastic medium) . aj(p) 
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saturates for large p a.t to^ — p/(0), yielding an enor- 
mous pile-up of extended (but non-propagating) states 
which causes the DOS to be concentrated at this fre- 
quency This is the glass analogue of a van-Hove 
singularity. At finite p, the density fluctuations of the 
x^'^ split this degeneracy, and yield the semicircular part 
of the spectrum at high frequency that can be recognized 
in the upper panel of Fig. 0. But the semicircular spec- 
trum misses the Debye part, and a better approximation 
is needed. So we substitute G in the last term of the 
r.h.s. by the resolvent of the continuum elastic medium 
Go{z,p) = {z — uj'^ (p))^^ . This is reasonable because the 
/^(g) factor makes low momenta dominate the integral, 
and due to translational invariance G{z,p) w Gi^{z,p) in 
this region We shall be looking at small so to a 
good approximation 



(27r)3 



1b 
p 



where the sound velocity is c 
positive constant. Then Eq. ( 
and can be solved to give 



(4) 



p/"(0)/2, while B is a 
is again quadratic in Q, 



pfiO) + B 
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[u^ - p/(0) + B + 1 



(5) 



We have two limiting cases. At high densities and low 
frequencies {pf{0) ^ w^, B, 2^pA), i.e. when the semi- 
circular part of the DOS does not reach low frequencies, 
the square root can be Taylor-expanded, and one gets 



2Trpc^ 



(6) 



which is precisely Debye's law. At small densities, on 
the other hand, the center of the semicircle (which is 
at = p/(0) — B) starts to be comparable to its ra- 
dius (cx ^/p), implying that the states in the semicircle 
hybridize with the sound waves. This is the source of 
the mechanical instability we find numerically in Fig. ^. 
Mathematically, the instability arises when G{0) develops 
an imaginary part. This can only come from the square 
root in Eq. (^), and it will happen for p < pc, with p^ 
fixed by the condition 2^ Apc + B ^ pc/(0). Now when 
/9 > Pc and uj <^ Lj* — 2-KC'^ pcAj (pc/^(0)), the square 
root in Eq. (^) behaves as 



^JW 



Pc 



(7) 



with D a positive constant. One distinguishes two 
regimes: 

• Lu*D{p — Pf.) <^ Lo ^ Lo*: the imaginary part of G is 
proportional to ^/uj, the DOS being g{uj) cx w^^^. 
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FIG. 2: Top: g{ijj)/uP from the numerical solution of the self- 
consitent equation of jl^. Bottom: Order parameter of the 
mechanical instability phase transition vs. density (points). 
Sohd line is a fit to the predicted behavior A(pc — p)^^'^ ■ 



• uj <ti U}*D{p — Pc) <C OJ*'. the imaginary 
part of Q is proportional to w, and ^(w) = 
wV(wVPc^£'(p- Pc)). The DOS is Debye-like 
with a prefactor arbitrarily larger than in Eq.(|^). 

One thus identify wbp with Lu*D{p — pc) that, on ap- 
proaching the instability, becomes arbitrarily small com- 
pared to any natural frequency scale. The mechanical 
instability is a phase transition, for which the order pa- 
rameter is — ImfJ(O). From Eq. (|^) we see that this or- 
der parameter behaves as (pc — p)^ , with (3 = 1/2, like 
in mean-field theories. In this sense, we can say that 
Eq. (^ is a mean-field theory for the model. That our 
qualitative analysis is correct can be checked looking at 
Fig. ^, where we show the DOS, and the behavior of 
—lmQ{0) with p, both obtained from the self-consistent 
G{q,z). Note also that the analytical solution compares 
reasonably well with the numerical solution (Fig. |l|) on 
a qualitative level (though it overestimates pc, and likely 
finds the wrong critical exponent). Since the behavior 
of the propagator at high momentum does not strongly 
affect the dispersion relation , we do not expect devi- 
ations from a linear dispersion relation at lubp ■ This can 
be checked either numerically (Fig. |^, bottom) or solv- 
ing numerically the self-consistency equations. This is in 
sharp contrast with lattice models for the BP |2^, p3| . 
Indeed, on the lattice the DOS is given by the integral of 
the dynamic structure factor, which is dominated by the 
Brillouin-peak. Thus the so-called lattice BP reflects a 
non smooth behavior of the Brillouin-peak. In fact, one 
can check (Fig. ||) that in these models the BP arise from 
a transverse van-Hove singularity that affects strongly 
the longitudinal dispersion relation. 

In summary, we have discussed an harmonic model for 
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FIG. 3: Dispersion relation of a heavily disordered (A = 1.0) 
FCC lattice model |2^] along the (1, 1, 1) direction, and (ro- 
tated) g{uj)/uj'^ , from the CPA approximation (a is the lattice 
spacing). The BP is at the end of the transverse branch. 



the high frequency dynamics of glasses, without trans- 
verse modes. The model has a mechanical instability 
transition controlled by the density. We have studied 
the model numerically, and analytically with a mean- 
field theory for the instability transition. The vibrational 
spectrum contains a BP which is the precursor of the 
transition. The BP in our model shares the main features 
of the experimental BP: it appears for frequencies in the 
linear part of the dispersion relation and it shifts towards 
arbitrarily low frequencies on approaching a mechani- 
cal instability. We also reproduces qualitatively the be- 
haviour of the silica BP when the density changes |l^ 
(the detailed theory of the temperature evolution of the 
silica BP should consider its negative thermal dilatation 
coefficient). The BP is built from the hybridization of 
sound waves with high frequency modes (extended but 
non propagating) that get softer upon approaching the 
instability. The analogous of our instability transition in 
nature is the topological phase transition Q that under- 
lies the dynamic crossover at the Mode Coupling temper- 
ature of real glasses |^ . The precise nature of the high- 
frequency modes that hybridize with the sound-waves is 
most likely material dependent and non-universal: they 
could be transverse 0, ^ |6), optical jl^, or even 
longitudinal modes as in our model. We believe however 
that the basic mechanism for the formation of the BP 
uncovered in our model is common to most (if not all) 
the structural glasses. Yet all real glasses do have trans- 
verse excitations, and one could ask about generic new 
features introduced by these modes. We are currently 
working to extend the present theory in this direction. 
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